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HELIUM RECOMBINATION LINES AS A PROBE OF ABUNDANCE AND 

TEMPERATURE PROBLEMS 



El resumen sera traducido al espanol por los editores. The paper presents a simplified formula to 
determine an electron temperature, T c (HeI), for planetary nebulae (PNe) using the He I A7281/A6678 line flux 
ratio. In our previous studies of T e (HeI) (Zhang et al. 2005), we used the He I line emission coefficients given 
by Benjamin et al. (1999). Here we examine the results of using more recent atomic data presented by Porter 
ct al. (2005). A good agreement is shown, suggesting that the effect of uncertainties of atomic data on the 
resultant T c (HcI) is negligible. We also present an analytical formula to derive electron temperature using 
the He I discontinuity at 3421 A. Our analysis shows that T e (HeI) values are significantly lower than electron 
temperatures deduced from the Balmer jump of H I recombination spectra, T (H I), and that inferred from the 
collisionally excited [O III] nebular-to-auroral forbidden line flux ratio, T c ([0 III])- In addition, T C (HI) covers a 
wider range of values than either T (HeI) or T c ([0 III]). This supports the two- abundance nebular model with 
hydrogen-deficient material embedded in diffuse gas of a "normal" chemical composition (i.e. ~ solar). 



The paper presents a simplified formula to determine an electron temperature, T c (HeI), for planetary nebulae 
(PNe) using the He I A7281/A6678 line flux ratio. In our previous studies of T c (HcI) (Zhang et al. 2005), 
we used the He I line emission coefficients given by Benjamin ct al. (1999). Here we examine the results of 
using more recent atomic data presented by Porter et al. (2005). A good agreement is shown, suggesting that 
the effect of uncertainties of atomic data on the resultant T (HcI) is negligible. We also present an analytical 
formula to derive electron temperature using the He I discontinuity at 3421 A. Our analysis shows that T c (Hc I) 
values arc significantly lower than electron temperatures deduced from the Balmer jump of H I recombination 
spectra, T c (H I) , and that inferred from the collisionally excited [O III] nebular-to-auroral forbidden line flux 
ratio, T c ([OlII]). In addition, T C (HI) covers a wider range of values than either T c (HeI) or T c ([OlII]). This 
supports the two-abundance nebular model with hydrogen-deficient material embedded in diffuse gas of a 
"normal" chemical composition (i.e. ~ solar). 
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RESUMEN 



ABSTRACT 



There are two long-standing problems in nebular 
astrophysics that are termed the "abundance prob- 
lem" and the "temperature problem" . The abun- 
dance problem refers to the findings that when the 
heavy element abundances are measured, particu- 
larly in planetary nebulae (PNe), the results de- 
rived from collisionally excited lines (CEL) are of- 
ten lower than those derived from optical recombi- 
nation lines (ORLs). One manifestation of the tem- 
perature problem is that electron temperatures de- 
duced from the collisionally excited [O III] nebular- 
to-auroral forbidden line ratio - hereafter T e ([0 III]) 
- are systematically higher than those determined 
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1. INTRODUCTION 



from the Balmer jump (BJ) of HI recombination 
spectrum - hereafter T e (HI) (see Liu 2004 for a re- 
cent review). Some recent attempts to solve these 
problems have implications that PNe may have com- 
plex physical conditions which are far from being 
understood. A correlation between the abundance 
discrepancies and the temperature discrepancies is 
found by Liu et al. (2001), suggesting that the two 
problems may have a common origin. Two possi- 
ble solutions are a) the presence of temperature and 
density variations within chemically homogeneous 
nebulae (Peimbert 1967; Vicgas & Clegg 1994) and 
b) a two-abundance nebular model with hydrogen- 
deficient material embedded in diffuse gas of a "nor- 
mal" chemical composition (i.e. ~ solar) (Liu et al. 
2000). To further investigate the problems, more 
methods to probe nebular physical conditions are ob- 
viously valuable. 
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In an earlier paper (Zhang et al. 2005, here- 
after Z05), we presented a method to use He I re- 
combination lines to measure electron temperatures 
of PNe - hereafter T c (He I) . We studied the effect 
that temperature and density variations inside neb- 
ulae have on the T e (HeI) value and found that for 
a chemically homogeneous nebula, the possible pres- 
ence of temperature and density variations causes 
T (He I)>T C (H I). In contrast, the two-component 
nebular model predicts T c (He I)<T C (H I). There- 
fore, a comparison between T e (He I) and T C (H I) pro- 
vides an opportunity to discriminate between the 
two paradigms. Applying the method to a sample 
of PNe, we found that T c (HeI) values are signifi- 
cantly lower than T e (Hl) values, in agreement with 
the prediction of the two-component (also called two- 
abundance) nebular model. 

This paper is an extension of Z05. In Section 2, 
we present an analysis formula to derive T c (Hc I). In 
Section 3, we give a statistical discussion of T e (He I), 
T (HI) and T c (OlII). 

2. METHOD 

In order to derive T e (HcI), we use analytic for- 
mulae for the emissivities of He I lines as a func- 
tion of electron temperature given by Benjamin et 
al. (1999). We suggested that the He I A7281/A6678 
intensity ratio serves as the best line ratio suitable 
for temperature determinations (Z05). A non- linear 
equation, I(7281)//(6678) = /(T e ), was used in Z05 
to determine T e (HeI). 

However, we find that when T e < 15000 K, there 
is a simple linear relationship between T e (HeI) and 
the He I A7281/A6678 intensity ratio, as shown in 
Fig. 1. With an electron density iV c = 10 4 cm -3 , a 
least-squares fit yields the following relation 

T c (Hc I) = 49300 x - 2150 K. (1) 

The relation is quite insensitive to electron density. 
We thus recommend equation 1 as a general formula 
to determine T c (HeI) for T c < 15000 K. Considering 
the uncertainties of electron density, the maximum 
systematical errors of T e (He I) derived by equation 1 
is 7%. In Fig. 2 we compare the T G (HeI) derived 
from equation 1 and those derived by Z05 using the 
non-linear equation mentioned above. An excellent 
agreement is indicated. 

Recently, based on improved atomic data, Porter 
et al. (2005) presented new calculations of He I line 
emissivities. They claimed that for the 32 He I emis- 
sion lines they considered, the average difference be- 
tween the new He I emissivities and those of Ben- 
jamin et al. (1999) is 4.6%. In Fig. 1 we overplot 
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Fig. 1. Electron temperature is plotted against the He I 
A7281/A6678 intensity ratio for an electron density of 
10 4 cm -3 . The solid line shows the prediction based on 
the atomic data given by Benjamin et al. (1999). The 
dashed line is a linear fit (see equation 1). The trian- 
gles represent the prediction based on the more recent 
calculation of Porter et al. (2005) . 

the values obtained according to Porter et al. (2005) 
(cf. their Table 1). A good agreement between the 
two sets of data is shown, implying that there is 
hardly any effect of the uncertainties of He I emis- 
sivities on the temperature determination from the 
He I A7281/A6678 intensity ratio. 

The determination of T c (HcI) is based on the 
assumption of Case B for the He I lines. An impor- 
tant consideration, therefore, is whether the possible 
departure from pure Case B to Case A recombina- 
tion for the He I singlet lines may affect tempera- 
ture determination. In Z05, we also used the He I 
A7281/A5876 intensity ratios to determine electron 
temperatures, which are found to be in good agree- 
ment with those derived from the He I A7281/A6678 
ratios. Given the argument that the deviation from 
Case B has a different effect on the He I A7281 and 
the He i A5876 lines, we suggest that the effect is 
negligible. 

Another method to measure T (HeI) is to use 
the He I discontinuity at 3421 A, produced by He + 
recombinations to the He I 2p 3 P° level. A relation 
between the He I discontinuity at 3421 A and electron 
temperature is derived as, 



T c (Hc I) = (2267 ±46) x 



A7(3421) 
7(3613) 



-3/2 



K, (2) 



where A/(3421)//(3613) = [/(A3421") - 
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Fig. 2. Comparison of the T e (Hel) derived from equa- 
tion 1 and those derived by Z05. The solid line is a y = x 
plot. 

7(A3421+)]//(He I A3613). The observation of 
the He I A3421 discontinuity, however, is very 
difficult due to its weakness. In this paper, we 
present only a fit relation between this discontinuity 
and electron temperature. The application of this 
method to PNe is beyond the scope of this paper. 
We note that to apply this method, high quality 
spectroscopic observations of this discontinuity are 
required. 

3. RESULTS AND DISCUSSION 

We calculate T c (HcI) for 48 PNe. Our results 
show that T c (HeI)<T o (HI)<r c ([0HI]). The aver- 
age values for the sample PNe are T c (HcI)= 6300 ± 
2100 K, T e (HI)= 10300 ±3100 K, and T e ([Om])= 
11900 ±2600 K. The result is exactly opposite to the 
predictions of the scenarios of temperature fluctua- 
tion and density inhomogencitics but in good agree- 
ment with the expectations of the two-abundance 
nebular model proposed by Liu et al. (2000). 

Fig. 3 shows the distribution of T c (HcI), T e (Hl) 
and T e ([OlII]). Inspection of Fig. 3 shows that ex- 
cept for a few extreme cases, T c ([0 III]) falls between 
9000 and 15 000 K, and T e (Hel) falls between 3000 
and 9000 K, whereas T e (HI) covers a wider range, 
varying from 5000 to 15000 K. In the scenario of the 
two-component nebular model, the H-dcficient com- 
ponent is extremely cold so that no CELs arise but 
ORLs are heavily enhanced. Therefore, T e (HeI) and 
T c ([OlII]) characterize the H-deficient cold gas and 
the diffuse hot gas, respectively Due to its hydrogen 
deficiency, the ultra-cold ionized gas has much less 















T e (He I) 




— 1 — hh — i — h 


-+ 


— i — 

r 


t— 




-4 


IT 






i i i i i i 

r e (H i) 

n. . . . n : 




-4 




r,([o in]) ' 

~h m n : 



5000 10000 15000 20000 

r, (K) 



Fig. 3. Histograms showing the numbers of PNe with a 
given T e (Hel) (upper panel), T c (Hl) (middle panel) and 
T c ([O III]) (lower panel). 

of a contribution to the H I recombination spectrum 
than the He I recombination spectrum. Accordingly, 
T C (HI) is a weighted average over both regions. As 
a result, T (HI) covers a wider range of values than 
T (He I) orT c ([0 III]). 

In previous studies of a large sample of PNe 
(Zhang et al. 2004), we have found that the dis- 
crepancies between T e (Hl) and T c ([OlII]) are anti- 
correlatcd with electron densities, i.e. high-density 
nebulae have the smallest temperature discrepan- 
cies. Robertson- Tessi & Garnctt (2005) found a 
negative correlation between the CEL/ORL abun- 
dance discrepancies and the nebula diameter and 
the Balmer surface brightness, suggesting that bright 
and compact PNe have small CEL/ORL abundance 
discrepancies. Given that the abundance discrep- 
ancies are positively correlated with the tempera- 
ture discrepancies (Liu et al. 2001), the two find- 
ings are mutually consistent. We thus infer that 
the abundance and temperature discrepancies are re- 
lated to nebular evolution; with expansion of nebu- 
lae, the contribution from the H-deficient material to 
recombination spectra becomes increasingly impor- 
tant. Here we compare the discrepancies between 
T c (HeI) and T (HI) against nebular densities and 
CEL/ORL abundance discrepancies. However, no 
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prominent correlation is seen, implying that some 
properties of the postulated H-deficient inclusions 
which dominate T c (HeI), such as its helium com- 
position, are very different for different PNe. To fur- 
ther understand the H-deficient inclusions, high S/N 
ratio, high spectral resolution spectroscopy of ORLs 
from heavy elements is required. 

The presence of H-deficient inclusions may have 
an important influence on the determination of the 
He/H abundance. Based on a large sample of PNe, 
Z05 estimated that the filling factor of H-dcficient 
components has a typical value of 1(U 4 , which may 
cause the He + /H + ionic abundance ratio to be over- 
estimated by a factor of ^ 1.25. The estimation, how- 
ever, is based on a very simple assumption, that the 
H-dcficicnt inclusions in all PNe have identical prop- 
erties to those given by Pequignot et al. (2003) in 
constructing two-abundance photoionization models 
of NGC 6153. Further spectroscopy of PNe should 
provide more constraints on the physical conditions 
of the H-deficient inclusions. 

YZ acknowledges the award of an Institute Fel- 
lowship from STScI. The work of YZ and XWL has 
been supported partially by Chinese NSFC Grant 
NO. 10325312. Support for RHR comes from the 
NASA Long-Term Space Astrophysics (LTSA) pro- 
gram. 

REFERENCES 

Benjamin, R. A., Skillman, E. D., & Smits, D. P. 1999, 
ApJ, 514, 307 

Liu, X.-W. 2004, in Planetary Nebulae beyond the Milky 
Way, eds. J. Walsh, L. Stanghellini and N. Douglas 
(ESO, Garching bei Miinchen), in press 

Liu, X.-W., Luo, S.-G., Barlow, M. J., Danziger, I. J., & 
Storey, P. J. 2001, MNRAS, 327, 141 

Liu, X.-W., Storey, P. J., Barlow, M. J., Danziger, I. J., 
Cohen, M.,& Bryce, M. 2000, MNRAS, 312, 585 

Peimbert, M. 1967, ApJ, 150, 825 

Pequignot, D., Liu, X.-W., Barlow, M. J., Storey, P. J., & 
Morisset, C. 2003, in IAU Symp. 209, Planetary Neb- 
ulae, eds. S. Kwok, M. Dopita, R. Sutherland (San 
Francisco: ASP), p.347 

Porter, R. L., Bauman, R. P., Ferland, G. J., & 
MacAdam, K. B. 2005, ApJ, 622, L73 

Robertson- Tessi, M., & Garnett, Donald R. 2005, ApJS, 
157, 371 

Viegas, S. M., & Clegg, E. S. 1994, MNRAS, 271, 993 
Zhang, Y., Liu, X.-W., Liu, Y., & Rubin, R. H. 2005, 

MNRAS, 358, 457 (Z05) 
Zhang, Y., Liu, X.-W., Wesson, R., Storey, P. J., Liu, Y., 

& Danziger, I. J. 2004, MNRAS, 351, 935 



HE I LINES 



5 



Y. Zhang: Department of Astronomy, Peking University, Beijing 100871, P. R. China; Space Telescope Science 

Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA (zhangy@stsci.edu). 
R. H. Rubin: NASA Ames Research Center, Moffett Field, CA 94035-1000, USA. (rubin@cygnus.arc.nasa.gov) 
X.-W. Liu: Department of Astronomy, Peking University, Beijing 100871, P. R. China, (liuxw@bac.pku.edu.cn) 



